EFFECTIVE INTEGRATION OF THE NONLINEAR 
VECTOR SCHRODINGER EQUATION 
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4. Differentials and Integrald 



J N ELGIN, V Z ENOLSKI, AND A R ITS 

Abstract. A comprehensive algebro-geometric integration of the 
two component Nonhnear Vector Schrodinger equation (Manakov 
system) is developed. The alhed spectral variety is a trigonal Rie- 
mann surface, which is described explicitly and the solutions of the 
equations are given in terms of ^-functions of the surface. The final 
formulae are effective in that sense that all entries like transcenden- 
tal constants in exponentials, winding vectors etc. are expressed in 
terms of prime- form of the curve and well algorithmized operations 
on them. That made the result available for direct calculations in 
applied problems implementing the Manakov system. The sim- 
plest solutions in Jacobian "iJ-functions are given as particular case 
of general formulae and discussed in details. 
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1. Introduction 

The Vector Nonlinear Schrodinger equation (VNSE) usefully models 
the propagation of a polarized optical beam along an optical fiber. The 
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vector nature of the dependent variable models the polarization state 
of the beam. It is intended in this article to derive and investigate a 
general class of periodic and quasi-periodic solutions of this equation. 
As the spectral curve for this system is trigonal - rather then hyper- 
elliptic as for the scalar case - existing formulae for these solutions 
are rather formal and not tractable for applications. Here, we use a 
method first devised by Krichever |Kri77j to effect an explicit integra- 
tion of the VNSE. This approach permits us to investigate some special 
cases where the formal solutions thus obtained reduced to simpler types 
expressible in terms of hyperelliptic or elliptic functions. 

In this paper we shall consider the integrable 2 dimensional focusing 
Vector Nonlinear Schrodinger equation (VNSE) 

It was proven by Manakov |Man74j that this system is completely in- 
tegrable and, in consequence, ()l.H1.2p are now known as the Manakov 
system. 

Manakov's method is based on the Lax representation 

(1.3) 0. = M0, 

(1.4) 0i = 50, 

where 

' -iz qi q2 

(1.5) M{z) =1 -qi iz 



-q2 i z 



and 
(1.6) 



X 



2^2 - |gip - |g2p 2iqiz-qix 2iq2Z-q2. 
B{z) = -i\ -2iqiz-qi^ -2^^ + |g^|2 q^q^ 

-2iq2Z-q2x qiq2 -2z^ + \q2\'^ 

where bar denotes complex conjugation. The Manakov system can be 
represented in the form 

(1.7) B,~Mt = [M,B]. 

The simplest solution Manakov's soliton has the form 

(1.8) Qsoii^^t) = 2r/sech(2r/(x + 4et)) 

X exp{-2i - 4i (^^ - r7^)t}c. 
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where c = (ci,C2)"^ is a unit vector, |cip + |c2p = 1, independent of 
both X and t, and ^ and 1] are real constants. 

Periodic and quasi-periodic solutions expressed in terms of explicit 
^-functional formulae have been quoted by several authors. The one 
component case, i.e. standard nonhnear Schrodinger equat ion was d e- 
veloped in |Its76j . |IK76j . and |Pre85j (see also monog raph |BBE+94j 'l. 



The multi-component case was studied in |Kri77t lAHHQOj . while the 
special case of reduction to a dynamical system with two degree of 
freedom was studied in |CEEKOO] . In recent years attention has been 
directed to modulation instabilities of the multi-component equation 
and searching for homoclinic orbits |FSW()()j . IFiVlMWOOj . |WF(K)j . Al- 
though we are not touching this interesting and important subject we 
believe that effective 6'-functional formulae could shed some new light 
on it. Indeed, we believe that they will be as useful for studying homo- 
clinic orbits of the Manakov model as the one component 6'-functional 
formulae are for studying the homoclinic orbits of the standard nonlin- 
ear Schrodinger equation (see Sections 4.4 and 4.5 of |BBE+94p . 

The article is organized as listed in Contents. The work is a mixture 
of analysis and computer algebra implementations using the Maple 
code described in |DvH01j . 

2. Zero-curvature representation 

Denote by ti = x,t2 = t, . . . , t„, . . . a set of "times" and introduce 
the set of 3 X 3 matrices Li{z) , L2{z) , (2), .. . satisfying the zero 

curvature representation, 

<9 d 

where the matrices £1 and £2 are chosen to satisfy the Lax representa- 
tion (jl.3p - More generally, ^n{z) is expanded as the n-th degree 
polynomial 

(2.2) i:„(;2) = (22fLo + (22)"-% + ... + (2^)L„_i + L„, n = 1, . . . , 
having the property 

(2.3) Lliz) = -L„(^), 

where dagger f denotes conjugate transpose: = c^. 
Here 



(2.1) —Ljiz) - —L,{z) = [L,{z),Lj{z)], 



-|i 0^ \ / 

^112;' ^'~\--q O2 
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= if 

where 

g = (51,92)^ 
while, for k > 2 introduce the following ansatz 

O A^ T _ ( OLk 01-1 ^ I r(0) 

In equation ()2.4|) ^1 denotes a 2 x 2-matrix and L^°^ denotes a constant 
matrix of the form 

/c^i 

(2.5) Lf = 4^ 4,3 

\ '^3,2 '^3,3 

with arbitrary entries Cp ^. In this article we set Lf = 0. 
The following theorem is valid 

Theorem 2.1. The entries to the matrix \2.4^ in the zero- curvature 
representation are defined as follows 

• the vectors /3„ and 7„ are given by the equations 
where D acts on a vector f{x) as 

X 

(2.6) Df{x) = ^f{x) + J {q{x')\f{x')}^dx'q{x), 

where {■,-}a denotes the matrix which is the anti-hermitian part 
of anticommutator, so that 

(2.7) {a\ b}A = [a^b - b^a) I2 + ba^ - ab\ 

Therefore, the flows are defined as 

(2.8) q^^-^^q=\{iDTq. 

• The (1, 1) element of the matrix is defined recursively as 
follows 

k k-2 

(2.9) ak+2 = - i 7Lj/3j - i ^1 «fc-i«j+2, 

j=0 j=0 



EFFECTIVE INTEGRATION OF THE VNSE 



5 



with 

i 

ao = --, ai = 0, a2 =iqq, 

while the associated right lower 2x2 minor Ak+2 is given re- 
cursively as 

k k-2 

(2.10) Ak+2 = i Yl 'yk-jf^J + i Yl 

j=0 j=0 

with 

i _ rp 

yio = -I2, ^1 = 0, A2 = -iqq. 
In each case, contributions from the second sum appear only for k > 2. 

Proof. The proof of these results follows from the substitution of ansatz 
()2.2|) . ()2.4p into ()2.1|) with i = 1 and j = n and solving the equation 
recursively. Also, one has to take into account that ^Tr(iL^) =0. □ 

3. The spectral curve 

The spectral curve is fixed by defining the stationary flow as fol- 
lows: let the system depend only on times ti, . . . ,tn-i- Then the zero 
curvature representation ()2.H) written for £1 and has the form 

(3.1) ^Ln{z)[L,{z),Ln{z)]. 

This relation suggests we consider the polynomial equation 

(3.2) fiz,w) = 0, f{z,w)=det{iin{z)-wh). 
We shall call the polynomial equation 

(3.3) X ■.= {{z,w)\f{z,w) = 0}. 

the spectral curve. Evidently coefficients of monomials z'^w'' of the 
polynomial f{z,w) are constants of motion. In what follows we shall 
consider the Riemann surface of the curve X, which we shall denote 
by the same letter. 

To proceed we recall that any rational function of its arguments, 
(f){z,w) is called a function on the curve f{z,w) = 0. The order of 
the function (f){z, w) on the curve X is the number of common zeros 
{zi,wi),...,{zn,wn) of equations f{z,w) = and (f){z,w) = 0. The 
curve is hyperelliptic if it admits a function of the second order, it is 
trigonal if it admits a function of third order etc. 

In the case considered, the spectral curve can be written in the ex- 
plicit form as 

X = {iz,w)\fiz,w) = 0}, 
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2n-l n-2 

(3.4) +(^ - ^(2.)") \,{2zf-'^~' + Y,f^n-2-,{2z) 

j=n j=0 



where 2n — 1 parameters Aj i = n, . . . ,n — 1 and fij, j = 0, . . . ,n — 
2 are constants of motion and can be taken arbitrary, but satisfying 
conditions given below in ()3.9|) . The coordinate z of the curve is a 
function of the third order and therefore the curve is trigonal. 

The parameters Xj of the curve X can be computed in terms of q as 
follows: 

^ n— 1 

k=j—n 
^ n-2 

(3.5) ~ 2 ^ ak+20:j-k-i, j = n, . . . ,2n - 1. 

k=j—n—l 

In particular, 

(3.6) A„ = aoa„+i, 

(3.7) A„+i = aoan+2 + ^llf^o + ^7o/3n- 

The structure of the second term in ()3.4j) has been obtained analyti- 
cally, including the stated expressions for \j. By contrast, information 
concerning the final term has been obtained using Maple, which gives 
the polynomial structure of degree n — 2 indicated. 

It follows from ()2.3p that the curve X admits the anti-involution 
property 

(3.8) a : X — > X, where a : {z,w) {z, —w) 

That implies in accordance with explicit formula for Aj 

Aj = Aj, i = n, . . . , 2n — 1, 
3.9 ' 

^^j = -^^j, j=0,...,n-2. 

Therefore we have 



(3.10) a o f{z,w) = f{z,-w) = -f{z,w), 

what means that the curve X has required anti-involution property. 
Let us clarify now the question on the genus g of the curve X. 

Lemma 3.1. Let Ln^ = and the curve X is given by the equation 
^3-4^ with parameters Xi, fij in general position. Then the genus of X 
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is given by the formula 

(3.11) g = 2n-3. 
Proof. Write equation ()3.4j) in the form 

f{z,w) = {w+^{2zr){w-^{2zrr 

(3.12) +{w- ^(2^)")P„_i(^) + P„_2(^) = 0, 

where Pn-i{z) and Pn-2{z) are polynomials of degrees n — 1 and n — 2 
correspondingly. The discriminant of ()3.12|) be of the form 

/ d 

Discriminant (X) = Resultant [f{z,w),-^f{z,w),w 

= 256iP„_2(^)^'" + 16P„_i(^)2^2" + 27P„_2(^)' 

+ 72iP„_2(^)Pn-i(^)^" + 16P„_i(^)2z2« + 4P„_i(^)3. 

The degree in z of the Discriminant (X) be 4n — 2 because coefficient 
of the leading power be 

(3.13) A2-4i/io^O, n = 2,3,... 

for Aj, /ij in general position. Moreover for general values of parameters 
the Discriminant (X) has no multiple roots and all zeros are simple 
branch points of the curve X. Beside of that we remark that the curve 
X has no branch points at infinities, ooi, 002, 003 Therefore the curve 
X has 4?2 — 2 simple branch points altogether, which we will denote ei, 
62, e4n_2 . The application of the Riemann-Hurwitz formula 

(3.14) g = ^-N+l, 

where B is total branch number, being equal in the case 4n — 2 and N 
is the number of sheets of the cover over Riemann sphere, which is 3 
in the case, completes the proof. □ 

We remark that our formula for genus ()3.1H) is addressed to the 
concrete curve which is fixed for our analysis. The inclusion of constant 
matrices L^^^ can increase the genus. The discrepancy of our formulae 
with results of [AH H90j and |Wri99j is due to the fact that in there an 
estimate of upper bound for genus was given for more general curve 
then our be. 

Introduce further the Riemann surface of the curve X. To do that 
we define local coordinate ^(P) of a point P = {x,y) E X in vicinity 



8 



J N ELGIN, V Z ENOLSKI, AND A R ITS 



of another point P = {z,w) E X as follows 

{z + ^ ii P = {z,w) is regular point, 
a + if P = (a, w{a)) is branch point, 
I if P = (oo, oo) is regular point at infinity. 

To comment this definition we remark that for general values of param- 
eters Aj and fii the curve has only simple branch points with ramifica- 
tion number one, what leads to the structure of the second line of the 
definition. The curve X has 3-sheeted structure with regular points at 
infinities, ooi, 002 and 003 where the coordinate of the curve behave as 
follows 



2"~^ i i 

w = — 2''^"^ + the first sheet, 

w = — h — (A„ + a/A^ — 4 i (Iq)^ + 0{^^) on the second sheet, 

w = — h -(A„ - J\l - 4i/io)^ + C'(^^) on the third sheet. 

We shall also assume that the branch points are all complex, form 
the conjugated pairs, i.e. 

(3.17) e2fc_i = e2fc, Im 62^-1 < 0, 

and Re 62^-1 < Re 62^+1. 

We are in position now to introduce a suitable homology basis on 
the Riemann surface of the curve ()H.4j) . A canonical basis of cycles 
ttj and bi respecting intersection property o = 0, bi o bj = 0, 
ttj o bj = —bi o aj = 6ij which also respect the involution property 

(3.18) a{aj) = -aj, 

a{bj) = bj - 2aj - ^a^. 

The homology basis for the case g = 3 is shown in figure 1; here, the 
solid, dashed and dash-dotted lines connecting points ei to 62 etc. are 
cuts connecting the first to second, second to third and third to first 
sheets respectively. See caption for further comments. The homology 
basis for higher genera can be plotted analogously. 
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Figure 1 . Basis of cycles of the curve X of genus 3. The 
sohd, dashed and dash-dotted hues denote paths on the 
first, second and third sheets respectively. Correspond- 
ingly the solid to dashed line, dashed to dot-dashed line, 
and dot-dashed to solid lines illustrate trajectories pass- 
ing through these cuts. The cuts are similarly encoded 
for clarity. 



12 



23 



31 



Figure 2. Contours passing from sheet 1 to sheet 2, 
from sheet 2 to sheet 3 and from sheet 3 to sheet 1 



4. Differentials and Integrals 

4.1. Holomorphic differentials and integrals. Let X be algebraic 
curve ()3.4p of genus g and let du{Q) = {dui{Q), . . . , dug{Q) Y be the 
set of canonical holomorphic differentials, which are given at n > 2 
explicitly as 
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1 z- 



(4.1) 



dz, J 



dw - 



duj{Q) 



.n-2. 



dz, j = n — 1, . . . ,2n — 3, 



where f{z, w) be the polynomial defining the curve ()3.4|1 . At n = 2 the 
curve is elliptic; this case is studied in detail in Section 6. 

If du = (j){z,w)dz is an Abelian differential, then the action of the 
involution a, a* say, is defined by the relation 

a*du = (f){z, —w)dz. 

For the differentials dui we have that 



a*duj 



(4.2) 

Introduce the matrix of a-periods 



-duj. 





( 


. \ 


A = 


/ 


duk 






) 



',k=l,...,g 

Prom the properties ()3.18|) and ()4.2j) of a* it follows that all elements 
of the matrix A are real. Normalized form of the above differential is 
introduced as 

9 

(4.3) dvj = y^^Cjidui, 

1=1 

where the matrix C = A~^. Evidently, 

(4.4) dv] = -a*dvj. 

Now introduce r-matrix, as a matrix of b-periods of normalized dif- 
ferentials, 









. \ 


(4.5) 


r = 


/ 


dvk 








J 



j,k=l,...,g 



Then using (jTTR|l and (jO]) we find 

Tj^= j) dvk = - <j> a*dvk = 

= -Tjk + 25jfc + y^^sik- 



dvk 
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In other words we have 

(4.6) T = -r + To, 

where all diagonal elements of tq are 2 and all off-diagonal elements 
are 1. 

4.2. Meromorphic differentials and integrals. Our construction 
is based on the existence of certain Abelian integrals fli{Q) and fl2{Q) 
of the second kind, and similar integrals h2{Q) and h^{Q) of the third 
kind. We shall define these integrals as follows: 

Definition 4.1. Define normalized Abelian integrals ^i{Q) and VL2{Q) 
of the second kind 

(4.7) I dn,=0, j dn2 = 0, j = l,...,g, 

where dfli{Q), dQ2{Q), are the second kind Abelian differentials in such 
a way that integrals Qi{Q), and fl2{Q), have poles at the infinities ooi, 
002, 003 in the vicinity of which the following expansions are valid 

(-iz + l(E, + E2) + 0{l) at OOl, 

(4.8) ni(Q) = 1(^2-^1) + 0(i) at 0^002 

[iz-l{E2-Ei) + 0{l) at g^oog; 

and 

{-2i^2_l(^^^^^) Q — , 

2iz^ -1{N2-Ni) at g^oo2, 

2iz'^ + 1{N2- Ni) at Q — > 003, 

where Ei, E2 and A^i, A^2 are certain constants. 

We now compute b-periods for the second kind differentials. Let P 
be a point in the vicinity of ooj and { be the local coordinate. Then 
the expansion of the vector of normalized holomorphic integrals reads 

p 

Ui + V^^^ + ^ _^(i)^3 ^ _ _ _ ^ i = 1, 2, 3, 



(4.10) J dv 



where Ui are constant vectors depending on the initial point Q and 



Ui^ J dv. 

Q 
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It follows from the Bilinear Riemann Relation written for the differ- 
entials dQi,dvj and dQ2,d'Vj, that 

(4.11) V= ivW- iV(2)_ iv(3), 

(4.12) W = AiW^^'^ - AiW^^'^ - AiW^''\ 
where the vectors V and W are defined as 



(4.13) V, = -^ <{dQ,{Q), W, = -^ (fdQ.iQ), 
Zm J ZiTX J 



J = l,•••,^• 



We shall refer below to the winding vectors V and W as the main 
winding vectors while the vectors V*-*^ and W'^'^\ i = 1,2,3 we shall 
call auxiliary winding vectors. 

It is easy to see that the differentials dQi and satisfy the same 
symmetry property (j4.4p as the differentials dvk- Hence, similar to the 
derivation of (j4.6|) . we arrive to the relations 

(4.14) V, = V„ Wj = Wj,...,j = l,...,g. 
Moreover, we claim that 

(4.15) E'k = -Ek, iVfc = -iV,,...,fc = l,2. 

To prove the symmetry properties ()4.15|) let us notice that the con- 
stants Ei^2 and A^i_2 can be determined via the asymptotic relations. 



(4.16) / dni = -2iz + Ei + 

IQ2 



1 



(4.17) j^^ d^li = -2iz + E2 + o(^ 



(4.18) / dfi2 = -4iz^ - A^i + ( - ) , z-^oo. 



z 

(4.19) j'^ dn2 = -4:iz^ - N2 + O (^^^ , z^oo, 

where the point Qj belongs to the j-th sheet of the Riemann surface 

T^iQl) = 7r(Q2) = TTiQs) = Z, 

and TT : X — > C is the canonical covering map. We assume that z is 
a large real positive number and that the contours of integration in 
the integrals J^^ and Jg^ do not intersect the basic cycles and do not 
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pass through the points ooj. We shall denote these contours [2 and [3, 
respectively. The involution a acts on the contours Ij as follows, 

a 

(4.20) (t([2) = l2 + Yl + ^^'^"1 + ^^'^"2 + k^^^as, 



u=l 

9 



(4.21) aik) = [3 + ^43)a, + m(3)a, + n(3)a3 + A;(3)«2, 

1^=1 

where aj denotes a positively oriented circle around the point oOj and 
^(2,3)^ ^(2,3)^ fc*-^'^-*, and nl^'^^ are integers. It should be also noticed, 
although we won't use it right now, that the integers m^'^'^^ and n*^^'^-' 
have the property^, 

(4.22) m^^') - ra(^) = odd number, j = 2, 3. 
From (jl^ and (ji^ we obtain that 



dQu = <t> dVtk = — (b a*dVLk = — (t dVti 



k 

Qj -^h -^h >^o"{fj) 



rffifc = - / rffifc, J = 2, 3, fc = l,2, 

and ()4.15|) follows in virtue of the asymptotic ()4.16p - ()4.19|) . 

Definition 4.2. Define normalized Abelian integrals h2{Q) and h^{Q) 
of the third kind 

(4.23) (f dh2 = 0, (f dh = 0, j = l,...,g. 

Joj Jaj 

The integral h2{Q) has logarithmic singularities only and only at 
infinities on the first and second sheets, 001^2 where it behaves locally 
as 



"'^In the 5 = 3 example featured in the figure 1 , the contour [2 starts at the point 
Q2 on the second sheet, goes below the branch points to the branch point ei, passes 
to the first sheet and goes below the branch points to the point Qi. The contour 
[3 starts at the point on the third sheet, goes below the branch points to the 
branch point 63, passes to the second sheet, goes to the branch point ei, passes to 
the first sheet and goes below the branch points to the point Qi. The specifications 
of the integers m(^'^\ n^^'^^ and /c^^'^^ in equations (|4.2Q(I and l|4.21(l are: 

nfUl, 42)=„f=0; m(2)=fc(2)=0, n^^^ = -I, 

nf^=nf^l, nf=0; m^^) = 0, n^^) = fc^^) = -1. 
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(4.24) h^iQ) 



In z — ln(52 + 0(1) at Q — ^002, 
— In 2 + ln52 + 0(1) at Q — > ooi. 



The integral h^{Q) has logarithmic singularities only and only at 
infinities on the first and third sheets, 001.3 where it behaves locally as 



(4.25) h^{Q) 



ln2; — ln(53 + 0(1) at Q — ^003, 
-In 2; + ln(53 + 0(1) at Q — > ooi. 

In fl4.24|) and ()4.25p 82 and ^3 are certain constants. 

Observe that the b-periods of integrals h2{Q), h-^iQ) are given as 

001 OOl 

(p dho = [ dvj, (D dh-i = [ dvj, 7 = 1, . . . , g. 

27ri4 J ^' 27ri^ ' J ' 

002 003 

Denote these periods as — r2,3, so that 

002,3 

(4.26) J dv = r2,3. 

001 

Taking into account that a{oOj) = ooj (cf. (j3.16|) ) and, once again, 
()4.4j) . we conclude that 

(4.27) = -r2,3 mod Z. 

The a - invariance of the leading terms of the asymptotic of the 
differentials dh2, dh^ at the points cxoj. (cf. (I4.24p and (j4.25|) ) implies, 
instead of (j4.4|) . the symmetry equations. 



(4.28) a*c//i2,3 = dh2,z. 

Similar to (j4.16j) - (j4.19p . the constant terms in the asymptotic (j4.24j) 
- (|4.25|) can be determined with the help of the relations 

(4.29) f dhk=(fdhk-2\nz + 2ln5k + o(-\ 
jQk J^k \^/ 

z +00, k = 2,3. 

Similar to the case of the integrals ^1,2, we can now use the symmetry 
properties (jOOl), dOH), and to see that 



(4.30) 52,3 = -S; 



2,3 
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dh): = (p dhk = (t a*dhk 



(p dhk = (p dhk + 271 im^''hes\ooidhk + 271 in^''hes\oo^dhii 



and ()4.30|1 follows in virtue of the asymptotic relation ()4.29|1 and the 
parity relation ()4.22l ) ( equation ()4.22|1 is now important - it is respon- 
sible for the minus sign in (j4.3(J|) ). 

In the next section we shall construct explicitly the integrals Qi 2{P) 
and /i2,3(-P) and compute the constants i?i,2, -^1,2, ^2,3 in terms of 9- 
functions of the curve X. 

4.3. ^-function and prime-form. The ^-function of the curve X 
with characteristic \e] 



^ El ... Eg 

is given by the formula 
(4.31) 

e[E] (v) = ^ exp { i 7r(n + £')^r(n + e') + 2 i 7r(n + eY{v + e)} . 



In this paper we are considering only half-integer characteristics, e'^, Ei 



or 2 for any k,l = 1, . . . ,g. Even characteristic [e] (e 



1) is 



nonsingular if ^[e](0) 7^ 0. Odd characteristic (e^*'^^^^' = —1) is nonsin- 



gular if the gradient V^[£](0) = (^^[e](«)|^^^, . . . , ^e[E]iv)lJ^ is 
non-zero. 

The canonical 6'-function is the 6'-function with zero characteristic 



(4.32) 



(v) = exp I i nnFrn + 2 i Tin^v^ 



The ^-function with a characteristic [e] possesses the periodicity 
property: 



(4.33) 



e[£](t; + efc) = exp{-2i7r4}^[£](t;), 

9[e\{v + Tfc) = exp{-i7rrfcfc -2i7!:vk - 2iv:Ek}9{E\{v) 
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where Ck = (0, . . . , 1, . . . , 0)^, = (ri^, . . . , Tgk)'^ and k = l,...,g. 
Using the relationship between r and r derived above (see (j4.6|) ) and 
definition of the 6'-function, we have that 



(4.34) e[e]{v)e-'''' 0[e]{v), 

and for the canonical 6 - function, 



(4.35) e{v) = e{v). 

The Schottky-Klein prime form |Bak95trFay73| is defined everywhere 
on X X X and is introduced by the formula 

e[e 

(4.36) £(P,g) = 




where P = {x,y) & X and Q = {z,w) & X are arbitrary points and 
^ [£]('«) is the 6'-function, with non-singular odd half- integer characteris- 
tic [e]. Concerning the characteristic [e] it is natural to suppose without 
loosing generality that the vector 

e = e + re' e (9) C Jac{X), 0(e) = 0, 

where {6) be ^^-divisor, is parametrized as 

(4.37) e = ^ dv- Kp, 

k=l ^ 

where Kp is vector of Riemann constants with base point P and points 
Pi, ... , Pg-i are different branch points of the curve X. 

The prime-form 8,{P,Q) vanishes only on the diagonal, P = Q, in 
the vicinity of which it is expanded in power series as 

(4.38) £(p, Q)^dap)dm = e(p) - aQ) + o(e(p) - e(Q)), 

where ^ (P) and ^ (Q) are local coordinates of the points P and Q around 

Po, e(Po) = 0. 

The prime-form ()4.36|) permits to construct symmetric second kind 
differential 2-differential which is called Bergmann kernel on X x X as 

(4.39) duj{P, Q) =—- In £(P, Q) dxdz 

oxoz 
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where [e] are non-singular odd half-integer characteristics. 

The differential duj{P,Q), where the coordinates P,Q are given as 
P = {x,y), Q = {z,w) has the properties: 

i) It is symmetric, du!{P,Q) = duj{Q,P). 

ii) It is holomorphic except on the diagonal set {P = Q) where it 
has a double pole. If the points P, Q are places in the vicinity of the 
point Pq and ^ is the local coordinate around Pq, ^(-Pq) = then the 
expansion of duj{P, Q) near Pq takes the form 

(4.40) 

^) = ( (^(p)_\(g))2 + ^^(^o) + o(e(p) - dmdm, 

where S{Po) holomorphic projective connection (explicit expression in 
terms of 5- functions is given e.g. in |Fay73| ) . 

iii) The a-periods taken in variable P or Q vanish, 

(4.41) ^duj{P,Q) = 0, t = l,...,g. 

Introduce following notations for the directional derivatives: 
dvf{v)=j2v,^f{v), 

k=l ^ 

etc., 

k=l 1=1 ' 

where V = (Vi, . . . , V^)'^, W = (Wi, . . . ,Wg)'^ are constant vectors 
and /(f) is a function of the vector argument v = {vi, . . . ,Vg)'^. The 
following theorem can be found in |Fay73| , [Jor92j concerning direc- 
tional derivatives along the ^-divisor. 

Theorem 4.1. For any nonzero vectors a = {ai, . . . , ag)'^ , b = (6i, . . . , bg)'^ G 
O and points Pi, ... , Pg-i on X the following identity holds 



(4.42) 




det [a\dv{Pi)\...\dv{Pg^i)] 
det [b\dv{Pi)\...\dv{Pg_i)]' 
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where the point e is given by 
(4.43) e = ^ j dv - Ki 



1 Pk 

9-1 -f 



fc=l 



and the matrices in \4-4'^ have been expressed by indicating each of g 
columns. 

Lemma 4.2. Let [e\ be non-singular odd half-integer characteristic of 
the curve X . Then 



(4.44) 



(4.45) 



[e] ijdvW^, 2^j = 1,2,3, 



9^«^[£](0)^0, 



1,2,3. 



det 



V^^^\dv{P,)\...\dv{P^_,) 



det 



dv{Pi 



Proof. The holomorphic differential dui = idz/f^ vanishes at ooj, i = 
1, 2, 3 to the order An — 8 = 2g — 2 and therefore {2g — 2)oOi is equivalent 
to the canonical class. Hence the vector of Riemann constants with the 
base point at ooj is a half-period |FK80j . The curve considered has only 
simple branch points (see the proof of Lemma 13. ip and therefore non- 
singular odd half-period corresponding to the characteristic [e] can be 
given by the formula ()4.43|1 with Pi, i = 1, . . . being branch points. 

First prove ()4.45|) for i = 1, j = 2. Suppose the opposite. Then the 
vanishing of the directional derivative dY(2)0{0) = will lead, according 
to ()4.42|) . to the vanishing of the determinant 

dv{oo2 

_ d^{002 

But there is no 002 among branch points Pi as that was shown ear- 
lier. The contradiction obtained proves the statement. Other cases are 
considered analogously. 

Prove ()4.44|1 . Suppose the opposite. Consider further the prime- 
form £(P, ooi) given by the formula (j4.36p . It is well defined because 
directional derivative 9^(i)6'(0) 7^ 0. According to principal property 
of the prime-form it vanishes only at P = ooi. Therefore the supposed 
vanishing of the ^-function should lead to vanishing of the directional 
derivative in the denominator. The contradiction obtained proves the 
statement. □ 

4.4. 6'-functional construction of meromorphic integrals. To con- 
struct the required second and third kind integrals Qi^2{Q) and /i2,3(Q) 
we first construct corresponding meromorphic differentials with the aid 
of prime-form introduced. 
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The normalized meromophic differential of the third kind, dh2 with 
the poles in ooi and 002 of the first order and residues ±1 in the poles 
is given as 



(4.46) 



dh2{P) = 4 In 



g(P,ooi) 

£(P,002)' 



Analogously the normalized meromophic differential of the third kind, 
dhs with the poles in ooi and 00 3 the first order and residues ±1 in the 
poles is given as 



(4.47) 



t{P, 003) 



We are in position now to give ^-functional representation for the 
second and third kind integrals, which permit us to compute 6 constants 
-E'1,2, Ni^2 and ^2,3 in terms of 6'-functions. 

Consider three quantities 




dv dz 



Q=OOi 



(4.48) 



k=l 



dvk 



p 

v= J dv 



P 



d^(i) In 9[s] I J dv\ , 



where i — 1, 2, 3 and the point Q has coordinates {z, w). 

Lemma 4.3. The quantities Qi\P), i = 1,2,3 are second kind Abelian 
integrals with unique pole of the first order at corresponding to the index 
i infinity oOi with a and b -periods, 



(4.49) 



dn?{P)=0, (t dn';>{P) = -2i7iVi^'' l^l,...,g 
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and following behaviour at the infinities ooi, 002, 003 on different sheets 
(4.50) 

(l + c?+OiO, (X2,, + O{0, 

n?iP) = Xi,2 + 0(0, nf\p) = {l + c? + 0(0, 

[Xi,3 + O(0, [X2,3 + O(0, 

'X3,l + O(0, 



1 , .(3) 



4'^ + 0(0, 



where 



(4.51) Xij = dyi.M9[e] f / h 2 ^ i = 1, 2, 3, 

^ ' ' "2 5^«^^[.](0) 9^«^^[5](0)' 

Proof. Relations ()4.49|) for periods follow immediately from definition 
and periodicity properties ()4.33|) of the ^-function, while expansions 
fl4.50|) result direct computation. The quantities ()4.5ip and ()4.52|) are 
well defined because of inequalities fl4.44|) and ()4.45|) . □ 

Similar statement is valid for three quantities 

(4.53) Q'i\P) = (29^(0 - 4„ lnd[e] I J dv 

where i = 1, 2, 3 and the point Q has coordinates {z, w). 

Lemma 4.4. The quantities ^l2\P), i = 1, 2, 3 are second kind Abelian 
integrals with unique pole of second order at corresponding to the index 
i infinity oOi with a and b -periods, 

(4.54) J dnf (P) = 0, j dnf (P) = -4 i TiWf^ l = l,...,g 
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and following behavior at the infinities on different sheets 
(4.55) 

r^ + c«+O(0, f>^2,i + O(0, 

^?iP) = \ 1^1,2 + 0(0, ^?\P) = U + + 0(0, 

[ri,3 + o(0, [^2,3 + 0(0, 

^3,1 + 0(0, 

fif (P)= <! 1^3,2 + 0(0, 

^ + cf +0(0, 



where 
(4.56) 




(4.57) 

.(0 



9^„^[£](0)V^ l_^v«,v«,vw^[^KO) ^ a^(,e[£](0) 



d^,.,d[e]{0) J 3 9^(.)^[£](0) a^w^[£](0)' 



Proof. The second term in ()4.53p does not contribute to periods - it is 
meromorphic function but the first term in ()4.53|) leads to the relations 
fl4.54|) . The expansions ()4.55p result direct computation. The quanti- 
ties ()4.56|) and ()4.57p are well defined because of inequalities ()4.44|) and 

(im. □ 



The normalized meromorphic differential of the second kind dfli{P) 
and dQ2{P) are then given as 

(4.58) ni{p) = -in[^\p) + inf\p) + inf\p) + Ci, 

(4.59) Q2{P) = ~2iQ'^^\P) + 2iQ'i\P) + 2iQ':^\P) + C2, 

where Ck, k = 1,2 are constants and integrals q[''\p) are given in 
flCT|l and (11321). 

The representations (j4.58l4.59p of the meromorphic differentials per- 
mits to compute the constants -Ei,2, ^1,2 and ^2,3 in terms of ^-functions 

Theorem 4.5. The following 9-functional expressions are valid for the 
constants -£^1,2,^1,2 one? 82.^: 

El = i (Xu + X21 + X31 - X32 - cf^ - cS'^ 

(4.60) ; 

E2 = i (Xi3 + X31 + X21 - X23 - - cl'^ 
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and 



(4.61) 



and 



iVi = -2i (Fi2 + K,i + F31 - - 4' 

iV2 = -2i (Fi3 + Kjl + >^21 - Y23 - 4 



(2) _ Jl) 
2 '-2 

(3) ^(1) 



(4.62) 52,3 = ^^^\/^^i.^sAem^J^^,.As]{0), 

where quantities Xi j,Yi j,^^ are defined in \4-51\l , \4-5t^ , \4-5<] ) and 

Proof. Consider first integral with the first order poles at infinities. 
Expand ()4.58|) at ooi, 002, 003 and compare with the asymptotic con- 
ditions ()4.8p to obtain equations 



1 „ 1 



- iXi2 + 1X32 + Ci + icf^ = ^E2 - i^i, 

- i Xi3 + i X23 + Ci + i cf ^ = ~E2 + i^i. 
We find ()4.6()j) and the following expression for the constant 

Ci = ^ (X12 + X13) - ^ (X32 + X23) + ^(cf ) + cf). 

Consider further the integral with the second order poles at infini- 
ties. Expand ()4.59j) at 001,002, C)03 and compare with the asymptotic 
conditions (j4.9j) . Solving linear equations as before we find (j4.6H) and 
the following expression for the constant 

C2 = ^ (yi2 + yi3) - ^ (F32 + F23) + ^(cf ) + cf ). 

Consider further the third kind integral 

(4.63) /,,(P)=in||^ + C.. 

t{P, 002) 

On the first sheet we have that, as P ^ ooi. 



e[e] I / dv 



ELl^^-.(002) 



hiP) = In r ' ' +C, + 0(0 



e[e] (^J dv^ \jTX=i 
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(4.64) =lne + ln E ^^^'^^ ^t^l^^ 

\ fc=i 



-ln^[£](-r2) + C, + O(0 
whilst on the second sheet, as P — 002, we have 



/i2(P) = -In^-ln, 



\ fe=i 



eW(o) 



\ k=l 

(4.65) +In^[£](r2) + C, + O(0. 

Comparison of expansions ()4.64|) and ()4.65p with the expansions 
leads to C, 



h = -f^ and 



(2) 



^[^](0). 



\ fc=i 



The expression for 62 in ()4.62p follows. The expression for ^3 in ()4.62|) 
is derived analogously. □ 

We emphasize that the constants described in the Theorem 14.51 are 
fundamental: expression for the constants "* coincide with accuracy 
to a trivial multiplier with values of the projective connection, S{P) 
(see fl4.40|) ) at infinities, S{oOi),i = 1,2,3. The quantities Xjj and 
Yij can be expressed in terms of multidimensional Kleinian ( and p- 
function whose classical and modern treatment, in the hyperelliptic 
case, can be found in |Bak95j and |BEL97j correspondingly. We also 
remark that analogous expressions for constants -Ei,2, ^1,2 in terms of 
^-functions and winding vectors for Thirring model which is associated 
with a hyperelliptic curve are obtained in [EGHOOI, see also (OHoHj. 

The important symmetry relations for the costants -Ei,2, ^1,2, and 
^2,3, which have been derived earlier (see ()4.15p and ()4.30|) ). can be also 
easily obtained from the 6'-functional formulae of Theorem 14 . 51 with the 
help of equation ()4.34p . 



5. Algebro-geometric solutions of the Manakov system 

We now summarize a list of basic objects which are related to the 
curve ()3.4p 

1. A homology basis of oriented cycles aj and as discussed in the 
Section 3. 

2. The differentials dvj introduced in the Section 4.1 are normalized. 
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3. The matrix of the b-periods of the trigonal curve X and the asso- 
ciated 6'-functions as defined by equations (j4.5p and (j4.3ip respectively. 

4. The Abelian integrals ^li{Q), ^2iQ), hi{Q) and /i2(Q), Q G X 
which are fixed by the conditions (Ol) . and KMi . 

5. An arbitrary divisor D with degree degD = g of general position, 



I.e. 



i=l 



where vr(P) is three sheeted covering 



TT : X 



TT OO 



[OOl, 002,003 J 



and Cj are the branch points of the curve X. 
The vector valued Baker- Akhiezer function 

*(g,x,t) = {MQ,x,t),tlj2iQ,x,t),^P3{Q,x,t)f 

is uniquely defined by two conditions. The first of these conditions 
describes the analytic structure of ^ on X/{ooi^2,3} 

I. ipiiQyXyt): i = 1;2,3 are meromorphic on X/{ooi,2,3}- Their 
divisor of poles coincides with D. 

The second condition describes the asymptotic behavior of ^'(Q, x, t) 
at OOl, 2,3 and shows that ^{Q,x,t) has essential singularities at 001,2,3 

II. As Q ^ C)Oi,2,3, the asymptotic behavior of ^{Q, x, t) is given by 
the equations. 



I +0{z~') 




at Q 



at Q 




exp(— i — 2 i z^t) 
exp( izx + 2i z^t) 



exp( i zx + 2i z^t) 



at Q 



003, 



where 62, G C are non-zero constants. Indeed, we shall take 62 and 
(^3 from ()4.24|) and ()4.25|) . respectively. 
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Then, ^{Q,x,t) = {ipi{Q,x,t),ilj2{Q,x,t),ip3{Q,x,t))'^ is uniquely 
determined by the conditions I. and II. and may be exphcitly con- 
structed by the formula 

(5.1) 

el J dv + r] e{D) 

MQ) = — -T exp{xn,{Q) + m^iQ) -Ex + m}, 

eljdv-D]e{T) 



ei J dv + T-r2]e{D-r2) 



I OOl 



^1 J dv-D]9iT) 

^OOl 



(5.2) X exp{xni{Q) + ^^]2(<5) - E'x + N't + /la}, 



e\ J dv + T-rs]e{D-r,) 

e(j dv-D]e{r) 



(5.3) X exp{xni(g) + tn2iQ) + E'x - N't + hs}, 
where 

(5.4) TVx + Wt-D 
and 

(5.5) ^ = £ / dv-K^,. 

Here, -FCooi is the vector of Riemann constants with the base point 
OOl. The constants E,E',N,N' are defined hj E = {Ei + ^2)/2, 
E' = -{Ei-E2)/2, and = (A^i + ^"2)72, A^' = -(iVi - iV2)/2, where 
Ei,E2, Ni and N2 are the basic constants from (j4.8j) and (j4.9p . 

The parameters appearing in the above expressions tlJi{Q,x,t) are 
defined in ()4.13|) . and vectors r2,3 are defined in ()4.26|) . 

The proof of formulae ()5.1l5.2l5.3p is based on the standard argu- 
ments of the theory of algebro-geometric integration: the Riemann 
theorem, which provides the condition I., the non-speciality of the di- 
visor D, which guarantee the uniqueness of the function ^'((5,a;,t), 
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and the periodicity properties ()4.33|1 of the 9 - function, which ensure 
that the equations (jS.lj) - (j5.3p define a single-valued (meromorphic) 
function on X. (For more details - see e.g. the similar proof for the 
usual, one-component NLS discussed in Ch. 4 of |BBE"'"94j .) 

We now fix some connected neighborhood U of the point 2; = oo on 
CP^ which has no branch points. Then, for each z & U, tt~^{z) contains 
exactly three points denoted by Qj E X, j = 1, 2, 3, so that Qj ooj 
when z — ^ 00. For z & U the matrix function 



(5.6) 



^'(z,x,t) = (*(gi,x,t),*(g2,a;,t),*(Q3,a;,t)) 



is now correctly defined to enable us to use the |RRE+94j version of 
Krichever's method |Kri77j to solve the VNSE. This will require the 
asymptotic form for \E'(z,a;, t), whose leading term is 



/ 



-2itz^} 



52 9(r2-r)e{D- 


r2) 


z e{D)e{r) 




g{— ix2— 2 i tz'^ 


} 






53 9{r-s-r)e(D- 





{—ixz—2 itz } 



xe 
xe 



{E2X-N2t} 



e{r2+T)d{D) 

e{r2-D)e{r) 

y^^{ixz+2itz^} 
^■^{-Eix+Nit} 



S2' 



^{ixz+2itz^} 



e(r2-r3+T)e{D-r2) 
e{r2-D)e{T) 

^Q{{E2-Ei)x-{N2-Ni)t] 



e(r3+r)e(£)) 

6(r3-D)e(T) ^ 

y^^{ixz+2itz'^} 

^Q{-E2X+N2t} 



e{r3-r2+T)9(D-r2) 
e{r3- D)e{T) 
^{ixz+2 i tz'^} 

XQ{-iE2-Ei)x+iN2-Ni)t} 



z_ {ixz+2\tz'^} 
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1 + 0(1) 



z 9(r2-D)e(r) 



5s e(r3+r)e(P) 



\ 



6>(r2-r)6>(£)-r2) 
z e{D)e(T) 
^^{Eix-Nit} 



^3 g(r3-r2+r)6>(D-r2) 

z e(r3-r»)0(r) 

Xe^-(^2-Si)a;+(Ar2-Afi)t} 



,53 6>(r3-r)6>(£)-r3) 

2 e(r>)6»{r) 



52 6>(r2-r3+r)e(D-r2) 
2 e(r2-D)6»(r) 

xe{(^2--E;i)x-(Af2-Afi)t} 



1 + 



X exp{ i zxJ + 2 i z^t J} 

where 

(5.7) J °" 



1 

z 
&2 



1 



2 



Theorem 5.1. Lei T) = Qi + . . . ,Qg be non-special divisor of degree 
g satisfying the reality condition, 

(5.8) D = D, D = J2 [ dv- Koo, 

■' OOl 

Then the solution of the Manakov system reads 



(5.9) 



r/ie constants Ei^2, ^1,2 and 62,3 are defined in ^.6(J{ ), ^ . 61 ) and (|^.6'ij| ), 



t/ie vectors r2,3 are defined in \4.2b\l , and the winding vectors V and 
W are given in (|-^.ii| ) - Ili4-^^ - 
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Following the methodology of |BBE"'"94j . we shall first prove two 
general lemmas. 

Lemma 5.2. Let \E'(2,x, t) &e 3 x 3 matrix function holomorphic in 
some neighborhood of infinity on the Riemann sphere smoothly depen- 
dent on X, t with the following asymptotic expansion at infinity 



I3 + ^mfc(x,t) 
fc=i 

24 



(5.10) X exp{izxJ + 2izHj}C{z), 

where J is defined in ^5. 7}) and Cx{z) = Ct{z) = Then assuming that 
\5.1U\) is differentiahle in x and t 

= M{z)+oiz-'), 
^t^-' = B{z)+o{z''), 

where z — > 00 and 

(5.11) M{z) = izJ + i[mi,J] = 



where 



and 



1 z q 

P izl2 



P = 2i ' , q = 2i ' 

"^1,31 / V "^1,13 



B{z) = 2iz^J + 2iz[mi,J] 

(5.12) +2i[m2, J] -2i[mi, J]mi. 

Proof. Direct calculations □ 

Lemma 5.3. Suppose that "^{z) satisfies the condition of the Lemma 
15. gl and equations, 

= M^, = Bm 

with M{z) and B[z) defined in i5.ll}) and i5.1^) . Then M(z) and 
B{z) are of the form presented in the Introduction, but p replaces q. 

Proof. lYiom \1>^ = it follows that 

(5.13) rrii^x = — i [^^2, J] + i [^1, J]mi. 
Given any 3x3 matrix A, we can represent it as 
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flu \ /a Ol2 CLl3 
A=\ 022 023 + ^21 \=Ad + Aoff. 
032 033 / \ OSI 

Note that [A, J](i = 0. To prove the Lemma we only need to check that 

/ • T • T 

Bo = 2i[m2,J]-2i[rm,J]mi= ( '.^ ^ 



Direct calculation shows that 



5o). = (-2iK,Ji]mO,= ( '^n^ T 



— i pq 

At the same time taking into the account 



iq 



□ 



We can now proceed with the proof of theorem 15.11 Consider the 
matrix \E'(2,x, t) defined in we claim that 

(A) ^(2;,a;, t) satisfies conditions of the Lemma [5.21 with 

/ 1 

C{z) =0^0 

(B) "^{zjXjt) satisfies conditions of the Lemma f5. 31 

Proof. (A) has already been established - see the asymptotic form of 
^'(z, X, t) presented above. Moreover, we have from this form the fol- 
lowing expressions for the relevant vectors q and p. 

(5.14) gi,2 = 2 i 62,3lP^^^^^^j^exp {-E.^^x + N.^^t} , 

^ 1^2,3 - > 



(5.15) P., = 2isJ-^:^:M^^0^^^^e.p{E,,. - iV,,i} . 

To prove (B) it is enough to show that ^x{Q) = M{z)'^{Q) and 
= B{z)'^{Q) for all Q e X. Consider the first equation. Put 



/(Q)=*.(g)-M(z)*(Q) 
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and note that if Q is in a neighborhood of the point ooj, then 

f(Q)^{Fiz))^, F{z) = <il,{z)-Miz)^iz), 7c{Q) = z, 
where {A)j denotes the jth column of a matrix A. We have at 2; ^ oo 

(5.16) F{z) = [^^.^~i-M]^ 

e 0(l)e-i 0(l)e-i 




Oiz~^)e 0{z)e-^ 0{\)e-^ 
0\z~^)e 0(l)e-i 0\z)e-^ 

where e = exp(— i — 2 i z^t\ 
It follows that 

0(1) exp{— i — 2 i z^a;} at Q ^ oo\, 
0(1) exp{ i 2;t + 2 i z^a;} at Q ^ oo2,'i- 



fiQ) 



Hence by the non-speciality of the divisor D w e conclude that f(Q) = 
(cf. |Kri77j : see also Corollary 2.26 |BBE+94j 1. which implies *^(Q) = 
M{z)'^{Q) and the validity of the first equation follows. 

The second equation can be proven by considering 

f{Q) = ^t{Q)-B{z)^{Q) 
and applying exactly the same arguments. 

□ 

As an immediate consequence we arrive to the following corollary. 

Corollary 5.4. The functions q{x,t) andp{x,t) defined in (|5.i^| ) and 
i5.1^) form a solution of the equations 

-^Pt+Pxx + '2Q'^PP = 0- 

Using conjugation properties dHHl), (|07|) . (lOUj) of the 

quantities V, W, r2,3, -Ei,2, A''i^2, ^2,3, the symmetry property ()4.H5j) 
of the ^-function, and taking into account condition ()5.8j) which we 
imposed on the divisor 2) we see that 

(5.17) pj{x,t) = ajqj{x,t), j = l,2, 

where 



(5.18) ai,2 = 



e{r2,3 - D) 



e{D) 
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Hence qi and q2 satisfy the evolution equations 



1 + + 2(ai|gi| + a2\q2\ )qi = 0, 



dt dx^ 
A trivial rescaling 



1 ^ + + 2(ai|'?i| + a2|g2| )g2 = 0. 



(5.19) qi ^ ^/alqi, q2 ^ v^^z 

complete the proof of the theorem. 

We emphasize that the solution ()5.9|) obtained is effective because 
computation of all parameters of solution, such as winding vectors, con- 
stants coming into exponentials was reduced to computation of holo- 
morphic differentials, their periods and ^-functions. These last com- 
putations are well algorithmized by Deconinck and van Hoeij in Maple 
for arbitrary curve, see also their paper |DvH01j . At the end of the 
paper we describe a computing procedure based on Maple software to 
compute algebro-geometric solutions to the Manakov system. 



6. Example: solution in elliptic functions 

In this section we shall show how the construction works in the sim- 
plest case of genus one. The spectral curve X reads 

(6.1) {w + 2iz^){w -2iz^f + {2\2Z + \^){w -2iz^) + ^q = Q, 
where the parameters A2, A3 and /io are 

A2 = "3 = i {llf^o + 7o /3i) = - i {(f^x - 

(6.2) A3 = a4 - i 7^/3o - i = qIQx + 
/^o = i \qixq2 - 52x91 

As commented earlier, the expression for /xq, and the fact that the last 
term in ()6.1|1 is a polynomial of degree zero is obtained from Maple. We 
also remark that parameters A2,3 are real whilst /xq be pure imaginary 
as that was stated in ()H.9j) 

Let us denote quantity (see (jH.lHj) ) A = A2 — 4i/io > 0. The dis- 
criminant of the curve has not multiple roots if and only if 

/ioA(27A2 - 64A^)[-27i/io(A + i/io)' + A^A^ + 54i/ioA - 270/io)] ^ 
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Figure 3. Basis of cycles of the elliptic curve X of genus 
1. The solid line, dashed line and dashed-dotted lines 
denote paths on the first, second and third sheets cor- 
respondingly. The cuts between the first and second, 
second and the third and third and the first sheets are 
denoted as correspondingly 

This curve is of genus 1 and admits the following behavior on the 
sheets at infinities 
(6.3) 

z — w — — — ^ — — + on the first sheet, 

2;=^, '"^ = ^ + ^(''^2 + V^)^ + 0{^^) on the second sheet, 

2;=^, w = ^ + j(A2- a/A)^ + 0(^^) on the third sheet. 
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The holomorphic differential on X be of tlie form 

(6-4) du = -. 

Denote A = f^du its a-period, tlien tlie normalized holomorphic dif- 
ferential dv = du/A. 

In the case considered we have r2 = —r^ = r, 

002 

1 



(6.5) ^ ~ ~A J 

OOl 

Indeed, because of expansions ()3.16|) the function on the curve w — 2iz'^ 
has second order pole at the first sheet and first order zeros on the 
second and third sheets. Then the Abel theorem says r2 + = 0. 

The auxiliary winding numbers ()4.13|) are computed with the aid of 
()3.16|) as follows 

yd) ^ yi2) ^ ^(3) _ -^/A + X2 

(g_g) 4A' 8A^/A ' 8Ay/A 

ly(l) = 1^(2) _ ^(2) 



Therefore we have for the main winding numbers 
(6.7) V = ^, W = 0. 

To perform further calculations it is convenient to transform the 
curve X to the form of standard Weierstrass cubic. Namely there 
exists a birational transformation T 

(6,8) 

T-^:X^X X = 4:i{w -2iz^), y = 32i z{w - 2i z^) + 8X2 

of the curve X to X = {x, y) 

(6.9) y"^ = - g2X - gs 
with parameters 

(6.10) (?2 = 64A3, g3 = -64A. 

The discriminant of X, coincides with multiplier 27 — 64 A3 of the 
expression for the discriminant X. In what follows we shall use Weier- 
strass functions of the curve X. 
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The transformation T maps canonical holomorphic differential du 
to dx/y. The point 2ujr appeared to be a zero of the Weierstrass p- 
f unction, 

(6.11) Pi^ur) = 

what follows from the fact that the map T is mapping ooi of the 
curve X to the oo of the curve X and 002,3 of the curve X to points 
(0, ± i ^/g3) of curve X. Remark that solution of 1)6.111) are given in 
|EZ82j in terms of Eisenstein series what will be of importance for 
developing computational algorithms. 
Direct computations gives^ 

(6.12) Ei = -E2 = E = -i(C(2cur) - 2r]r), N^^2 = 0. 

The equalities W = 0, Ni^2 = are in accordance to our analysis 
which leads to the statement that the time t = t2 is stationary for the 
flow yielded the Manakov matrix B (see 11.7)1 . 

The solution has the form 

a (r)-2iA ^3(V-x±r-D) 

Ai,2 = ± — ;= 1= exp < larg ' 



'4V2A VA Mr) I ""KM^r-D) 

where 6{z) = i^siz) is canonical ^-function of the genus one curve ()6.ip . 
D is arbitrary constant satisfying 6{D) 7^ 0. Substitution of the elliptic 
solution ()6.13p to the expressions for levels of the integrals of motion 
(16.2)) leads to equivalences and permit moreover to compute the link 
between periods A and 2u as^ 

(6.14) A = -2u. 

It is remarkable that the simplest solution to the Manakov system, 
i.e. solution of genus one independent in time and as the result it 
soliton limit enable to yield Manakov soliton p. 8)1 . We can guess fll.8)) 
can be obtain as the result of degeneration of the genus three curve. 

7. Summary: computational algorithm 

This section is addressed to a reader who wants to know the com- 
puting algorithm without going through the arguments of the paper. 



^We shall give this computation for completeness in the Appendix 
^see Appendix for details 
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The procedure to compute algebro-geometric solutions to the Vector 
Nonhnear Schrodinger equation 

can be formulate as follows. 

• Fix positive integer n G {2, 3, . . .}. 

• Fix a polynomial in two variables (algebraic curve) 

f{z,w) = {w + ^{2zr){w - ^{2zrf + {w- ^(2z)")P„_i(z) + iP„_2(^) = 0, 

with arbitrary real polynomials Pn-2 and Pn-i of degrees n — 
2 and n — 1 correspondingly. Compute its genus g by using 
|DvHnij and Maple. For polynomials Pn-2 and Pn~i in general 
position g = 2n — 3. 

• Compute the vector of holomorphic differentials du{Q) = {dui{Q), . . . , dug{Q))' 



1 Z' 

dujiQ) = -Q—, -dz, j = 0, . . . , n - 3, 

duM = a u j = n-2,...,2n-3. 

• Compute vector of normalized holomorphic differentials 
dv{Q) = A~^du{Q) 

• Compute auxiliary winding vectors V^'^\W^^\ Z^^\ i = 1,2,3 
from expansions 

p 

dv =0(1) + + W(^)e^ + Z«e^ + ..., i = 1,2,3, 

Set for main winding vectors 
V = i V^^) - i _ iy(3)^ = 4iVF(i) -4iVF(2) - AiW^'^K 

• Compute vectors r2,3 

0O2,3 

^2,3 = / dv. 



OOl 
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• Compute 6 constants 



■ (2,3) . (1) 



N^^2 = 4i (9^(1) - 9^(2,3)) ln^^[£](r2,3) - 4i9^(3.2)ln [e[£](r3,2)^^[£](r2,3 - 



+ 2 i 9^(3,2) ^(2,3) In 



0[e]ir3,2 



0[£]{r2,3 - r3^2 



(2,3) . ^(1) 



ICn — IC. 



2 ) 



where [e] be non-singular odd characteristic, d\r, dw and dz 
are directional derivatives, 

Sv-E"'-^. Sl^w = ttv^W,^ etc. 

fc = l fc=l 1=1 ' 

and also 6 constants, at i = 1, 2, 3, 

(,) ^ _ / g^,)g[5](0) \ ' l ^v(0,v(0,vW^N(0) g^(.)g[5](0) 

• The algebro-geometric solution is of the form 

, , BiVx + Wt- D + r23) , ^ , 

gi,2(x,t) = 2iAi,2 — g (^3; + vFt - ' ''''^ i--^i.2X + Ari,2t} , 



M,2 = '52,3exp < i ar 



(D) 



Mr2,3-D) 

where D is arbitrary vector satisfying condition 6{D) 7^ 0. 



8. Appendix 

Write in addition to ()6.6p the third auxiliary winding numbers 
(8.1) 



16 AAi 

The constants cl''\ I = 1,2 and k = 1, 2, 3 are 

(8.2) 4^)==0, _ A3(A2^±v/A) 
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and 



(i)_ 1 ^'i"(0) 



^2 48v42 ^[{0) ' 
(2,3) ^ i/ioAi(i/io-A^±yA) 1 (A2± VA)^^r(O) 

Compute first Ei — E2. We have 



2i V« (ln^9i(r))' - + v^{3)) (ln^9i(2r))' 



El - ^2 = i (X12 - X13) - i (X32 - X23) + i (ci'^ - c\' 

i Aq 



2VA 



Apply the duphcation formula 

C(2.) = 2C(.) + 1^, 
2 ^'(z) 

which in the case considered reads 

2Aq 

a^ur) = C(2cur) - 

to obtain 

i A3 /cj 1 



(8.3) Ei-E2 = V{AujC{2ujr) - 8r]ujr) ^ , . ■ 

Analogously obtain 

El + E2 = 2 i (X2 + X31) - i (X32 + X23) - i (cf ^ + cf ^ 
= i (^(2) - r(2))(2 (lni9i(r))' - (ln^?i(2r))') - 

i A2A3 f LU ^ 1 



i A2A3 
2A 



To complete computation of -^1,2 we must find relation between period 
A of the curve X and period 2uj of the curve X. We shall do that by 
substituting ()6.13|1 to the ()(i.2j) which should lead to equivalence. We 
have 

qixq2 - q2xqi = -4(5253exp[(£'i + E2)x] 



X 



n/ A2 (^i(^^' ~ r)^9'i(Vx' + r) - i9i(Vx' + r)^9'i(Vx - r)) 

Vl \/ X j 
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E2 — El 



■jj^^iiVx' ~r)^i{Vx' + r) 



where x' = x — ^{1 + r). Substituting instead of derivatives i!^'i{z) 
Weierstrass ^-functions we transform the expression in the curly brack- 
ets to 

[2.Va2Vx u + 2ur) 

- 2ujVC{2Vxuj - 2ujr) - 8Vujr]r + E2 - Ei] 



a{2ujrYexp 



2 

r rj 



p{2ujVx') 



X 



Au}VC{2ujr) - W'qujr + E2-E1 



p{2uVx') 



where we apphed addition formula for the Weierstrass a and (^-functions 
and took into account (|6.1H) . Because this quantity should be a con- 
stant with respect to x we set 

(8.4) El- E2 = 4:UjV{C{2ujr) -2rir), Ei + E2 = 

what in combination with ()8.3p gives relation ()6.14|) between a-periods 
of the curve X and X. Therefore the derivation of expressions for Ei_2 
given in ()6.12p is completed. But let us continue the computation of 
the integral level. We have now 

i3i(rY 

151x^2 - g2xgir = 256 '"""^ 



But 



6161 



i/io Moy 



^ di{<df 
and (73 = —64 A, 



256A4A-(?i(r)2 

what completes the derivation of equivalence. 

It remains to show that = 0. Develop expression for A^^i. It is 
the sum of 3 terms which are 



^2 ^1 



riuj( 



-3A| + 10 i /io + AaVA) i /ioAK i /io - + Aa^A) 



8A^A 



A^{\2-VAy 



2(1^(1) - 1^(2) + iy(3))(in^9,(r))' - 2W'^^\lnMMy = ^^kp- 
2((VW)2 + (V(2))2 + (F(3))2)(in^?,(r))" - 2(V^(='))2(lni9i(2r))" 
r/cu(-3A2 + 10i/io + A2VA) uj^XliX2 + VAf 



4A2A2 
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Taking into the account equality (jfi.llj) and substituting then flfj.l4|l to 
the sum we obtain necessary equahty. The equahty A'"2 = is derived 
in analogous way. 

It could be also of interest to check by direct substitution that so- 
lution ()6.13p satisfies to the Manakov system. To check that we first 
compute 

(8.5) 

2{\q,{x,t)\' + \q2{x,t)\') = 8{6l + 6l)^ ^-^^ 

i^[{0)^3{Vx' + r)^3{Vx' -r) a{2ujVx + 2uor)a{2uoVx - 2uor) 
~ §^{VxYdi{ry 2ct2(2 * uoV)a^{2uor) 

= -^p(2a;yx), 

where we used again in the Weierstrass addition theorem Fur- 
ther the first derivative 



d_ 

dx 



qi{x,t) 



qi{x, t) 2VujC{2ujr) - Wrju -E + 



Vuj{p'{2ujV)-i^) 



p{2ujVx) 

The first 3 terms in brackets vanish because of expression for E given 
in ()().12|1 . Therefore 



(8.6) 



dx"^ 



qi{x,t) 



qi{x,t) 



V p{2ujVx) 



+ 



d f Vuj{p'{2ujV)-i^) \ 
dx \ p{2ujVx) ) 



-qi{x,t)p{2ujVx). 



Combining ()8.5|) and ()8.6|) we obtain the equiality Vahdity of 

the ()1.2|) is proved analogously. 
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